Combinatorial influence of environmental temperature, obesity and cholesterol on SARS-CoV-2 infectivity

The continuing evolution of SARS-CoV-2 variants not only causes a long-term global health concerns but also encounters the vaccine/drug effectiveness. The degree of virus infectivity and its clinical outcomes often depend on various biological parameters (e.g., age, genetic factors, diabetes, obesity and other ailments) of an individual along with multiple environmental factors (e.g., air temperature, humidity, seasons). Thus, despite the extensive search for and use of several vaccine/drug candidates, the combinative influence of these various extrinsic and intrinsic risk factors involved in the SARS-CoV-2 virus infectivity has yet to be explored. Previous studies have reported that environment temperature is negatively associated with virus infectivity for SARS-CoV-2. This study elaborates on our previous findings, investigating the link between environmental temperature and other metabolic parameters, such as average total cholesterol and obesity, with the increase in COVID-19 cases. Statistical analysis conducted on a per country basis not only supports the existence of a significant negative correlation between environmental temperature and SARS-CoV-2 infections but also found a strong positive correlation between COVID-19 cases and these metabolic parameters. In addition, a multiphase growth curve model (GCM) was built to predict the contribution of these covariates in SARS-CoV-2 infectivity. These findings, for first time, support the idea that there might be a combinatorial impact of environmental temperature, average total cholesterol, and obesity in the inflation of the SARS-CoV-2 infectivity.


Results
Geographical locations of countries with highest levels of average total cholesterol, BMI, obesity, and COVID-19 cases in their population. To investigate whether the geographical location has any impact on SARS-CoV-2 infection, we first identified the top 75 countries for COVID-19 infection (Fig. 1A). As of December 2020, there are 75 countries with at least 14,000 COVID-19 cases, providing a relatively large number for comparative analysis between various parameters. Moreover, 40 out of these 75 countries belong to a relatively colder temperature zone i.e., 15 °C or below (average annual temperature) i.e., above 23.5° N latitude and towards the north pole, which suggests that areas with lower environmental temperature are more prone to SARS-CoV-2 infection. We simultaneously identified the top 75 countries for relatively high average total cholesterol, high-BMI, and prevalence of obesity in their population in the world map in order to associate the influence of these particular metabolic parameters with SARS-CoV-2 infections.
In our recent publication, we investigated whether the novel coronavirus may exhibit increased infectivity in regions with lower environmental temperatures 5 . Additionally, our laboratory had previously reported a negative correlation between environmental temperature and serum cholesterol content of the associated population, based on geographic locations 9 . Furthermore, it was found that the locations with the highest COVID-19 caseloads overlapped with those with the greatest number of metabolic parameters occupying the same locations ( Fig. 1B-D). Therefore, these data indicated that these metabolic variables, i.e. high-cholesterol, high-BMI, and high-obesity, might additionally affect viral infectivity.
Univariate analysis to study the association between COVID-19 cases and metabolic parameters. As stated earlier, our previous publication identified a negative correlation between environmental temperature and the number of novel coronavirus cases 5 . In this earlier publication, the COVID-19 case data www.nature.com/scientificreports/ was limited to the months of March and April of 2020 5 . Our findings were consistent with other studies [10][11][12] , detecting significantly lower COVID-19 transmissibility in warmer temperatures. However, some conflicting studies did not find any significant correlation between temperature with COVID-19 cases 13,14 . Generally, such studies were carried out in using only a single country/region. Thus, for this current study, we have expanded our analysis to include data for the total COVID-19 cases per million populations from March, 2020 to July, 2021. Univariate analysis using the Spearman correlation method found a strong negative correlation between MAET and the number of COVID-19 cases in the winter months in comparison to the other months for this period of time, thereby implying that the winter months may promote viral infectivity. A similar trend is also seen when using Pearson correlation analysis. Furthermore, we examined the association between the previously-discussed metabolic parameters and COVID-19 cases. The results of the univariate analysis between these metabolic parameters (average total cholesterol, BMI, and obesity levels) with COVID-19 cases per million are shown in Table 1. In agreement with our previous work, we again found a negative correlation between environmental temperature and COVID-19 cases 5 . Thus, to examine the interdependence of cholesterol levels and total COVID-19 cases, we conducted a univariate analysis between average total cholesterol and the month-wise COVID-19 cases per million by both Pearson and Spearman's correlation methods. A strong positive correlation was observed between cholesterol and total COVID-19 cases per million of a country (Table 1). Finally, we performed a univariate analysis between COVID-19 case number and BMI of the population of a country and found a positive correlation between these parameters (Table 1). Therefore, the metabolic parameters, i.e., highcholesterol and obesity, could contribute to novel coronavirus infectivity. However, univariate data analysis alone is not sufficient to address the combinatorial influence of these parameters on virus infectivity.
Analysis of the distribution of total COVID-19 cases per million from March 2020 to July 2021. To better understand the potential connection between the various factors and the number of COVID-19 infections, further exploratory data analysis was conducted. Figure 2, contains a set of box plots showing the distribution of COVID-19 cases per million over the specified months. These plots show that the mean and median both increase over the period, the mean is larger than the median indicating that the distribution of cases per million is positively skewed for the sample of considered countries. Additionally, the dot points in the initial months indicate that the data points from some countries were outliers. The difference between median and mean increased gradually even in the year 2021, clearly indicating that countries with extreme values were still skewing statistical measures. In addition to this, a simple trajectory plot for the monthly COVID-19 data from March, 2020 to July, 2021 ( Fig. 3) showed that in the early months of the year 2020, there is instant growth in COVID-19 cases per million in most of the countries included, in the middle of months of 2020 a stability is seen in all countries and finally, in the last few months, there is again rapid growth in most countries. It can be seen that the rising trend in cases among the countries is non-linear and the growth rate varies from country to country over the defined period.
Modelling approaches for examining the relationship between environmental temperature, average total cholesterol, and BMI and COVID-19 cases of a country. We used a latent growth curve model (GCM) approach to capture the growth of COVID-19 cases for considered time periods. Furthermore, the multiphase GCM elected to study the intra-individual and inter-individual changes in trajectories for considered countries. In the absence and presence of various covariates, we assigned the ranks to different fitted multiphase models as per their performances. It was observed that model-8 had the lowest AIC and BIC values and large enough TLI value (Table 2), and so it was used for modelling the COVID-19 case number of all countries' data. The estimates with corresponding standard errors obtained from model-8 are reported in Table 3. From (Table 3), it can be observed that the average baseline value for the outcome variable C[t] n is 65.4469 ( g 0n ). The average growth amount in the second, third, and fourth phases (as defined in the "Methods") are −0.0004 ( g 1n ), −1.5834 ( g 2n ) and 3.9195 ( g 3n ), with significant variances. The estimated parameters of vectors and A 3 [t] ( α 6 ,α 7 , α 8 , α 9 , α 10 , α 11 , α 13 , α 14 , α 15 , α 17 , α 18 ) can be used to calculate the value of outcome variable C[t] n at any month using the estimated values of intercept and slope. The variance terms explain the extent to which countries differ in intra-individual changes, and covariance terms are responsible for the interindividual differences among countries. Variances differing from zero indicate that these countries differ in their initial levels and slopes of cases per million. For the predictors, it can be inferred that one unit change in BMI is associated with −0.0193 ( γ 0 ) unit change in slope ( g 0n ), 0.0002 ( γ 1 ), −0.3019 ( γ 2 ) and 0.5621 γ 3 unit changes in intercepts ( g 1n , g 2n and g 3n ). Similarly, one unit change in CH is associated with 0.5393 ( η 0) unit change in slope ( g 0n ), −0.0010 ( η 1) , 1.0465 ( η 2) and −1.1209 ( η 3) unit changes in slopes ( g 1n , g 2n and g 3n ). For TVC, one unit change in temperature in March is associated with −0.0816 ( τ 3 ) unit change in C[t] n . Other coefficients τ 4 , τ 5 , … , τ 19 can be interpreted similarly. Here, it can be observed that TICs are responsible for the changes in intercepts and slopes. They have an indirect effect on C[t] n mediated by latent variables g 0n , g 1n , g 2n and g 3n , while TVCs have direct impact on outcome variable C[t] n . A structure plot has also been constructed for the outcome variable, cases per million, alongside the other predictors (Fig. 4). In the structure plot, C [3], C [2], …, C [17] represent the outcome variables while T [3], T [2], …, T [17] represent the MAET. The latent variables g 1 , g 2 , and g 3 represent the change in the outcome variable in the second, third, and fourth phases respectively. The structure plot depicts the association of temperature with the corresponding month's outcome variable with coefficients ( τ t ). Furthermore, the coefficients of BMI ( γ t ) and cholesterol ( η t) illustrate an indirect effect on the outcome variable through latent variables g 0n , g 1n , g 2n and g 3n . The dotted lines in the plot show the fixed parameters, while the bold line corresponds to the estimated coefficients from the model. www.nature.com/scientificreports/

Discussion
The novel coronavirus has spread to pandemic proportions, posing a major threat to the human population with an infection rate that has been found to increase exponentially. This outbreak of COVID-19 has led to millions of deaths worldwide. With the lockdown lifted in many parts of the world, infection rates have been found to increase rapidly. In fact, virus infectivity has rapidly peaked as recently as December 2021. There are multiple environmental and biological factors that could contribute to the rapid spread of disease, such as environmental temperature and metabolic parameters. Studies have reported an association between the infection and transmission of viruses with air temperature and humidity. This has been shown, for example, for the influenza virus 15 . Reports in the literature suggest that temperature might be an important factor accounting for the transmission of other coronaviruses, like SARS-CoV 16 and MERS-CoV 17 because of (i) increased virus half-life at lower temperatures, (ii) greater stability in nasal passages when the epithelial surface is cold, and (iii) greater stability in lower humidity as compared to intermediate humidity [15][16][17] . Considering these facts, we conducted a study to establish a link between environmental temperature and COVID-19 cases. We observed that countries with www.nature.com/scientificreports/   www.nature.com/scientificreports/ more COVID-19 cases were mostly located north of the latitude of Wuhan, China where the pandemic started in December, 2019. Thus, we performed a detailed country-wise statistical analysis which established a significant negative correlation between COVID-19 cases and MAET of a country 5 . However, this initial finding was limited to the COVID-19 data from March and April, 2020. While there are a few studies whose findings are broadly congruous with our own regarding temperature and SARS-CoV-2 caseload 18 , others have found no correlation between temperature and infection rate 13 . In this study, we have validated the relationship between MAET and COVID -19 cases per million from March, 2020 to July, 2021 of a country. Univariate analysis by both Pearson and Spearman's methods indicated a negative correlation between temperature and COVID-19 cases ( Table 1). The statistical analysis also identified a stronger negative correlation for the winter months (November to March) signifying that the warmer months faced fewer SARS-CoV-2 infections compared to the colder months. The box plot (Fig. 2) and the simple trajectory curves (Fig. 3) also demonstrated a high prevalence of COVID-19 infection in the later months of the year as compared to the initial months. To further reconfirm that lower temperatures influence infection rate, we inspected the geographical locations of the countries with moa greater number of Here, C represents the outcome variable; T represents the monthly average environmental temperature (MAET); CH represents cholesterol; g 0 , g 1 , g 2 and g 3 represents the latent variables. The dotted lines show the fixed coefficients and the bold lines show the estimated coefficients. www.nature.com/scientificreports/ infections and observed that most of these countries were located above 23.5 o N latitude and/or towards the poles, further suggesting that cold temperatures may affect the SARS-CoV-2 transmission (Fig. 1A). Nevertheless, we did not limit our study exclusively to environmental temperature since there may be multiple reasons for this increase in the novel coronavirus infections. More severe viral infection is observed in those patients who are already suffering from the other pre-existing health complications 19 . Therefore, we refined our research by considering additional metabolic parameters like high-cholesterol, BMI, and obesity, in conjunction with the environmental temperature of a country, in influencing the SARS-CoV-2 caseload. The role of cholesterol in virulence of other respiratory viruses, like influenza, has been well established in several studies. For instance, cholesterol, which maintains membrane structure, is critical for viral stability and virulence 20 . Studies have documented that patients with prior high cholesterol levels are more prone to viral infectivity, eventually leading to severe disease outcomes 21 . Cholesterol-enriched lipid rafts might accommodate the aggregation of ACE2 receptors on the cell membrane, thus enhancing the binding of the S-protein of SARS-CoV-2 to the host cell surface 22 . Another study has shown that individuals with an apolipoprotein (apo) E4 genotype have an increased risk of severe COVID-19 infection. Though increased cholesterol levels promote ACE2 and furin trafficking inside host cells 21 , cholesterol plasma levels are found to be decreased in patients post-infection. In brief, high cholesterol present in host cell membrane, virus particles, and human blood may augment the virus entry processing in the host cells 21,23 . Our data depicted that geographic locations above 23.5 o N latitude and towards the poles had a higher prevalence of high average total cholesterol levels, often overlapping the areas with the highest COVID-19 cases (Fig. 1B). Furthermore, the univariate analysis also showed a significant positive correlation between average total cholesterol and COVID-19 total cases per million (Table 1), thereby suggesting that higher cholesterol levels may enhance the infection rate of SARS-CoV-2.
Obesity is a critical health condition which is a consequence of modern sedentary lifestyle. In addition to its other health implications, there is an association between obesity and critical viral infections 24 . Reports have shown that obesity can contribute to the progression of viral infections such as in the case of Hepatitis C infection 25 . Several studies have reported that overweight patients need respiratory support and have increased admission to intensive care units (ICUs) compared to patients with normal weight, even at a younger age 26 . A cohort study showed that obesity is an important factor in disease severity of SARS-CoV-2, having the highest impact on patients with a BMI ≥ 35 27 . Moreover, in vitro experiments have shown that ACE2 and TMPRSS2, two essential entry components for SARS-CoV-2 infection, are highly upregulated in the lung epithelial cells of obese patients 28 . Therefore, we explored the relationship between BMI and obesity and COVID-19 caseload. Upon identifying the geographical locations of countries with a higher prevalence of BMI and obesity, significant overlap was observed with those with high numbers of COVID-19 cases (Fig. 1C,D), just as with cholesterol levels. Additionally, statistical analysis identified a positive correlation between the BMI with the total COVID-19 cases per million (Table 1).
Based on these preliminary findings from the univariate analysis data, we conclude that these metabolic parameters i.e., average total cholesterol and BMI, influence the infectivity of the SARS-CoV-2 virus. To authenticate these findings, different statistical approaches were used, and we attempted to model the COVID-19 cases/million trajectory using a latent growth curve model in the presence of time-variant and invariant factors. Multiphase GCM was used to investigate the role of metabolic parameters on the escalation of COVID-19 cases. We evaluated different covariates such as temperature, average total cholesterol, and BMI and tried to fit it in the multiphase models individually with the COVID-19 cases per million and examined the AIC, BIC and TLI values ( Table 2). The model having the lowest AIC and BIC values and greater TLC values was considered to be the best fit model fulfilling all the statistical criteria, and based on this, the models were ranked accordingly. Moreover, when these metabolic parameters were evaluated with environmental temperature, it was found to have a greater impact on the infection rate. In order to determine the combinatorial effect of all these factors i.e., environmental temperature, average total cholesterol, and BMI with the COVID-19 cases per million, we incorporated these parameters all together in the multiphase model and determined the AIC, BIC, and TLI values; it was observed that the AIC and BIC values are lowest and TLI is highest for this model in comparison to the other competing models. Therefore, it could be concluded that Model-8 outperformed the other models. The estimates used were obtained from this model. The structure plot for the data including all countries is shown in Fig. 4. Altogether, these findings indicate that patients with higher cholesterol, BMI, and obesity may be more prone to infection, particularly in the winter months. Thus, obese individuals with high average total cholesterol may be at additional risk for getting SARS-CoV-2 infection if they are further exposed to cold environment. This is the first attempt to model COVID-19 cases/million trajectory using the latent growth curve model in presence of time-variant and invariant factors. In fact, such a growth curve modeling approach could be utilized to track and predict the spread SARS-COV-2 infection over the time in presence of the considered factors. This can be helpful to design the policies against the COVID-19.
Although our study indicates a negative correlation between temperature and the number of COVID-19 cases, the ability of this virus to infect might also depend on age, sex, and ethnicity, the prevalence of different diseases in the population, different social distancing practices, and uses of various preventive medicines. Also, our findings are based on the effect of atmospheric temperature on COVID-19 cases; how indoor temperature might affect infection rate is yet to be considered. Our study has considered a rather holistic approach to understanding the role of temperature in infection rates of the virus and takes into account the fluctuations observed in a single country. Moreover, with the emergence of SARS-CoV-2 mutant strains, it is somewhat difficult at this stage to speculate the role of temperature, obesity, and cholesterol on the infection rates of these mutant strains. Although prior SARS-CoV-2 infection protects most individuals from reinfection for at least five months 29 , the first case of COVID-19 reinfection after recovery has been identified in a female from Japan 30 , after which reinfection became a true threat. A recent study has shown that people above 65 years of age have relatively low protection against reinfection by COVID-19 31  www.nature.com/scientificreports/ associated with lower antibody titre 32,33 . All these findings suggest that various metabolic factors not only enhance the infectivity rate, but also provoke reinfection. However, how metabolic parameters like obesity and cholesterol levels affect the incidence of reinfection still remains unclear. Having said so, the pattern of infection by the virus may differ in the near future due to our growing knowledge of treatment and a much-improved understanding of the SARS-CoV-2 virus infectivity and its associated complications. Furthermore, this study suggests that individuals with metabolic disorders such as high-cholesterol and obesity could be more susceptible to SARS-CoV-2 infection in the winter months, especially while living in a colder environment. ACE-2 expression in host cells and average total cholesterol levels may be increased in response to exposure to a cold environment and in the winter months 9,34,35 . Elevated ACE-2 levels were observed in multiple metabolic disease conditions such as obesity, diabetes, and higher LDL cholesterol 36 . Additionally, obesity may not only increase ACE-2 and TRMPSS2, but also cellular cholesterol levels by increasing SREBP1 37,38 . Interestingly, a role for low temperature has been suggested in stabilizing the RBD-ACE2 interface and triggering "open" conformations of the COVID-19 spike protein, thus enhancing viral infectivity at cold temperatures 39 . Several studies have also reported multiple roles for cholesterol in enhancing the susceptibility to SARS-CoV-2 infection. Cholesterol-rich microdomains can provide an effective platform for interaction between ACE2 and Spike S-protein 7 . Tang et. al reported the role of cholesterol in increasing the density of ACE2 receptors on host cell membranes 23 . Reports using super-resolution imaging have also observed increased SARS-CoV-2 entry in cells with high blood serum cholesterol levels 40 . Furthermore, studies have shown that obesity is a critical factor in COVID-19 severity. Increased ACE-2 expression in lung tissue is seen in patients that are obese, implicating excess adipose tissue in enhancing the spread of the virus 41 . Thus, a colder environment and obesity both increase ACE-2 and host membrane cholesterol, which favour viral entry processing resulting in increase of virus infectivity.
However, other intrinsic factors like hypertension, cardiovascular diseases, renal diseases, and cancer, as well as extrinsic factors like relative humidity and indoor temperature, have not been included in these analyses. Specifically, indoor temperature may contribute to virus infectivity. Actually, cold temperatures and low relative humidity (RH) adversely increase the half-lives of the virus 15 . Aerosolized SARS-CoV-2 has the potential to stay infectious for about 16 h at optimum indoor meteorological conditions 42 . Indoor locations have a relative humidity < 40% which indicates higher chances of airborne SARS-CoV-2 transmission. Moreover, SARS-CoV-2 infection occurs in cool, dry, air-conditioned indoor environments 43,44 and during cold weather, people mostly stay indoors, which further potentiates transmission 45 . Presently, the USA is setting up policies which suggest maintaining indoor temperatures between 20 and 24 °C and RH around 20-60% 46 . Thus, dry weather generated to maintain indoor temperature (20 °C) in winter months may further potentiate virus infectivity since the virus may persist for longer times in a relatively lower humid atmosphere 47 .
In brief as a conclusive remark, such multiphase growth curve models may be used to depict the contribution of various covariates with COVID-19 cases. Moreover, the individuals with metabolic disorders such as high cholesterol and obesity may have additional risk for this virus infectivity especially in winter months or while living in colder environments. Thus, this study further recommends that a nationwide policy is to be framed in order to combat COVID-19 pandemic and its clinical outcomes for taking care of vulnerable individuals with such metabolic diseases of a cold country. However, further study is required to know whether the infectivity rate of mutant variants of SARS-CoV-2 depends on these factors.

Methods
Data collection. COVID-19 case data was collected for every month from March, 2020 to July, 2021 from2021 from ourworldindata.org/coronavirus 48 . We considered the total cases per million as the variable of interest for the COVID-19 case dataset for each month. The environmental temperature data for different countries was collected from climatestotravel.com, as described in our recent publication 5 . The latest data for cholesterol and obesity were acquired from the WHO website and BMI from ghdx.healthdata.org. These datasets are publically available. All data collection was performed in accordance with the relevant guidelines and regulations.
Geographical distribution of countries. The geographic location of different countries has been marked in the world map with the help of world map marker (worldmapmaker.com), explained previously 9 . This online server is freely available. The top 75 countries for SARS-CoV-2 infections have been marked in the world map. The top 75 countries for highest average total cholesterol, BMI, and obesity levels in their population have also been marked on the world map.
Univariate analysis. Univariate analysis was performed using the Spearman and Pearson methods to determine whether to accept the null hypothesis between the two variables, as described previously 5,9,[49][50][51] . Statistical analysis was considered to be significant for p < 0.05 . Univariate analysis was performed using the software GraphPad Prism 6.0.
Statistical models. For this work, we applied the multiphase latent GCM to study the variable of interest (cases per million) over the defined time period by incorporating time invariant covariates (TICs) and time varying covariates (TVCs). Here, total cases per million over the months was taken as variable of interest while average total cholesterol and BMI were considered as TICs and temperature as TVC. Further, TICs and TVC were examined to see whether the changes in cases per million over months may be explained by the presence of these predictors. R software was used for this analysis. For the modelling purposes, "lavaan" package in the CRAN directory was utilized. www.nature.com/scientificreports/ In model fitting, a latent GCM approach was used for COVID-19 cases per million in the presence of different associated covariates: temperature, cholesterol and BMI. The latent GCM captures the growth trajectory of the variable of interest over the defined time period. TICs and TVCs were considered to see whether the changes in cases per million over months may be explained by the presence of these predictors. Here, cholesterol and BMI were treated as TICs, and temperature as TVC. The general structure for the GCM in presence of TIC and TVC can be given as below: where, C[t] n denotes the cases per million for different countries for any particular month t and vectors A 0 [t] , A 2 [t] , …, A k [t] are responsible for capturing the intra-individual trajectories of countries over the months. The latent variables g 0n , g 1n , …, g kn represent the slopes and intercepts based on which countries differ interindividual and τ t , γ k , and η k represent the coefficients associated with the covariates temperature, BMI and cholesterol respectively. Also, ϑ k 's are means of slopes and intercepts; ξ kn 's are residual deviations from latent slope and intercepts with some variances and finally ǫ[t] n denotes time specific residual with variance σ 2 ǫn . The TICs associate directly with the slope and intercepts and have an indirect effect on repeated measures of outcome variables while TVC poses a direct impact on outcome variables. Based on the above equation, several models like linear, quadratic, exponential, and latent GCM have been tested and we found the multiphase GCM to be the most suitable model. Multiphase growth curve model. The multiphase GCM is based on two or more regression lines which allows modelling of multiple processes that are responsible for the intra-individual changes over the time periods 52 . As observed from the trajectory plot, the COVID-19 case growth rates were different in the considered months. Therefore, it was more suitable to fit multiple splines than just one line for the defined time slot to capture the COVID-19 cases for different countries. Based on this proposition, we divided the total COVID-19 cases per million over the months into four phases for modelling purposes. Four vectors A 0 , A 1 , A 2 and A 3 respectively, were used to model the initial or baseline, second, third, and fourth phase. Since, we had the data for seventeen months from March, 2020 to July, 2021, we defined these four phases as Phase I: (March, 2020-May, 2020), Phase II: (June, 2020-December, 2020), Phase III: (January, 2021-April, 2021) and Phase IV: (May, 2021 -July, 2021). The vectors for the chosen multiphase model can be defined as: The parameter α i 's was estimated from the data. Moreover, the inter-individual differences in the extent to which each process contributes can be governed by the random vectors g 0 , g 1 , g 2 and g 3 . The variances among g 0 , g 1 , g 2 and g 3 represent the extent to which countries differ on each of the particular aspects of intra-individual change invoked by the corresponding vectors, A 0 , A 1 , A 2 and A 3 while the covariances among these variables were interpreted as the inter-individual differences in one aspect with respect to another aspect.
Some well-known fitting criteria such as AIC, BIC and TLI were used for choosing the best model. The lower values of AIC and BIC with sufficiently large value of TLI directs to choose a better model. So, a model having lowest values of AIC and BIC with large enough TLI was preferred to others.